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In this supplementary materials section, we provide (i) additional information on the numerical 
simulations of the main work, (ii) the derivation of all prefactors in the analytical theory, and (iii) 
unpublished experiments of the contact stiffness of a polymer pressed against a rough substrate. 

INTRODUCTION 

In this supplementary paper we present some details which could no be included in the PRL because of space 
limitations. We first show the surface roughness power spectra used in the simulations. Next we derive the prefactor 
in the power law relation between the contact stiffness and the load. Finally we describe a new measurement of the 
contact stiffness where a silicon rubber block is pressed against an asphalt road surface. 

NUMERICAL DETAILS 

Fig. Q] shows a surface roughness power spectrum as used in the simulations. The solid lines indicate the mean 
values for the spectrum, while the dots reflect one particular realization. Fluctuations of the height h(r) in real space 
are not only the consequence of variations in the absolute value of their complex Fourier transforms h(q) but also due 
to the random phases. 

DERIVATION OF PREFACTORS 

Consider a randomly rough surface with a roll-off as indicated in Fig. [2l The power spectrum 

C(q) = C for q < q < q r 




FIG. 1: (Color online) Power spectra for two surfaces without (a) and with (b) a roll-off at large wavelength as generated by 
a Fourier filtering algorithm. The solid lines show the prescribed power spectrum C(q) and the dots the actual realization. 
Panel (b) indicates the wavevectors of the long- wavelength roll-off q r = 7r/A r and the short-wavelength cutoff qi = 7r/Ai. For 
q < qo — tv/L where L is the linear system size the surfaces have zero power. The noise at low q is due to the fact that order 
q 2 Fourier components contribute to the power-spectrum of a realization of a surface. 



where go = if/L, where L is the linear size of the studied system. We also write q r = ir/A r where A r is the roll-off 
wavelength. The surface mean square roughness amplitude 
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FIG. 2: The surface roughness power spectra as a function of the wavevector (log-log scale) for a self-afEne fractal surface with 
a roll- off. 

We first calculate the elastic energy stored in the deformation field associated with the Hertz mesoscale asperity 
contact region. The mesoscale asperity has the radius of curvature R and the radius of the (apparent) contact region 
between the mesoscale asperity and the flat countersurface is denoted by ro. The mean summit asperity curvature is 
given by [l[ R = /3y / 2Ko where kq is the root-mean-square curvature of the surface 



*l = \J d 2 q q^Ciq) = IT jT dq q 5 C(q) 



When roughness occurs on many length scales so that qi/qo » 1 Nayak has shown that (3 — (8/37T) 1 / 2 . Including 
only roughness components with wavevector q < tt/tq gives the mean summit curvature 1/R of the mesoscale asperity: 
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We define R — q r R, h° ms — q r h° ms and ro = q r ro and assume 7r/ro > q r . Thus 
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From Hertz theory 



Defining 



gives 



R = ar, 



2-H 




3FR 

1e* 



F = 



Fgl 

E* 



= *FR = ^Ff, 



2-H 



Thus 



or 



3a 



11 = " ' t) 



(2-iT)/(l+iT) 



i? 



(2-^0/(1+^0 



_p(2-H)/(l+H) 



V 3/(1 + H)^(2-H)/(1+H) 



The elastic energy stored in the Hertz mesoscale deformation field: 



where S = q r S. We define 
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Next we calculate the elastic deformation energy stored in the vicinity of the microasperity contact regions in the 
Hertz mesoscale contact region 0, Q : 
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where U\ = q r u\. We have u\ = 7/i rms where 7 « 0.4 and 
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Using the definition for a to eliminate /i^ ms we can write 
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The total elastic energy 
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The total stiffness k = KAq is given by 
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The stiffness per unit area if = fc/Ao = k/L 2 : 
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In Fig. [3] we show l/0o, l/0i and 1/0 = l/0o + l/0i as a function of the Hurst exponent if. It is interesting to note 
that as if — > 0, then 1/0q — > while l/0i remains finite, i.e., for the fractal dimension Df = 3 the stiffness is entirely 
determined by the short-wavelength roughness in the macroasperity contact region. Note also that since q r f=a ir/ L, 
where L is the linear size of the system, the stiffness scales as k ~ q r H ^ 1+H ^ n*> L H /( l + H ) with the size of the system. 
This is in contrast to the region where the p ~ exp(— u/uq) relation holds where the interfacial contact stiffness is 
independent of the size L of the system. Note also that the stiffness scales with the rms roughness as {h® ms ) 
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FIG. 3: The quantities 1/0 O , 1/01 and 1/0 = 1/0 O + 1/01 are denned in the text. 



while in the region where the p ~ exp(— u/uq) relation holds the stiffness is proportional to (h° ma )~' s ~. For the Hurst 
exponent H ks 0.8, which is typical in practical applications, 8 « 1, which appears to be in good agreement with the 
prefactor found by Pohrt and Popov in their numerical simulation study [H . The treatment presented above can be 
generalized to obtain the distribution of stiffness values (at least approximately) by calculating the distribution P(R) 
of summit curvature radius R (which is easy to do). 

It is interesting to determine the critical force F c such that for F < F c one needs to use the finite size region 
expression for the stiffness while for F > F c the Persson expression is valid. When the relation p ~ exp(— u/uq) is 
valid the stiffness 
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The critical force F c is determined by the condition that k given by (1) and (2) coincide. This gives 
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Note that typically (for H w 0.8) (Q'y)^ 1+H )/ H « 0.2. The prediction (3) for the switching between the finite size 
region and the region where the stiffness is proportional to the loading force is in good agreement with our simulation 
results. To show this let us first write 



The surfaces we have studied in numerical simulations have the rms slope 0.1. To relate this to q r h° ms which enters 
in (4) we use that 
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In the present case the rms slope is 0.1 and qo/qi = 1/4096 and H = 0.7 so that q r h° ms ~ 5.5 x 10~ 3 for q r /qo = 
L/X r = 1, and q r h° ms « 1.3 x 10~ 2 for q r /qo — 8. Using this from (4) we get p c /E* w 8 x 10~ 5 for q r /qo — 1 and 
p c /E* «4x 10~ 6 for q r /qo = 8, which is in good agreement with Fig. 1 in our paper. For the surface with H = 0.3 
from (5) we obtain (for a surface with the rms slope 0.1) q r h° ms nearly 100 times smaller than for H = 0.7, which will 
shift the cross-over force F c , between the two stiffness regions, with a similar factor to lower values, again in good 
agreement with the numerical studies. The results presented above differ from the conclusion of Pohrt and Popov 
who state that the power relation observed for small applied forces is valid for all applied forces 0, In particular, 
in Ref. [5( Pohrt et al state: "It is the authors strong belief that the proportionality found by Persson appears only in 
the first case described. Whenever the surfaces are truly fractal with no cut-off wavelength, a power law applies." The 
present study shows that this statement is incorrect and Fig. 1 in our letter clearly shows that the contact stiffness 
cannot be described by a power law for all applied forces as this would correspond to a straight line on our log-log 
scale. 

As an example consider applications to syringes, where the relation between the squeezing pressure p and the 
(average) interfacial separation u (which determine the contact stiffness) is very important for the fluid leakage at the 
rubber-stopper barrel interface. Consider the contact region between a rib of the rubber stopper and the barrel. The 
width of the contact region (of order w p=> 1 mm) defines the cut off wavevector q r = ir/w ps 3000 m _1 . The Hurst 
exponent H ps 0.9 and the rms roughness amplitude (including the roughness components with wavevector q > q r ) 
is h° ms ps 3 /Lim. The elastic modulus of the rubber stopper is typically E ps 3 MPa. Using these parameters we get 
from (4) : p c « 1 kPa, which is negligible compared to the pressure in the contact region between the rib of the rubber 
stopper and the barrel, which is typically of order ~ 1 MPa. 



EXPERIMENTS 



The relation (1) as well as the above mentioned finite-size effect region has also been observed in experiments. 
In these experiments a rectangular block of silicon rubber (a nearly perfect elastic material even at large strain) is 
squeezed against hard, randomly rough surfaces. In this case no plastic deformation will occur, and the compression 
of the rectangular rubber block, {p/E')d (see below), which will contribute to the displacement s of the upper surface 
of the block, can be accurately taken into account. Such measurements were performed in Ref. Q, and were found 
to be in good agreement with the theory (these tests involved no fitting parameters as the surface roughness power 
spectrum, and the elastic properties of the rubber block, were obtained in separate experiments). Here we show the 
result for the contact stiffness K = —dp/du (not presented in Ref. |6j) of one additional such measurement. 

The experiment was performed for a silicon rubber block (cylinder shape with diameter D = 3 cm and height 
d = 1 cm) squeezed against a road asphalt surface with the rms roughness amplitude 0.63 mm and the roll-off 
wavelength Al f=a 0.3 cm as inferred from the surface roughness power spectrum. The squeeze-force is applied via 
a flat steel plate and no-slip of the rubber could be observed against the steel surface or the asphalt surface. We 
measured the displacement s of the upper surface of the block as a function of the applied normal load. Note that 

s = (uc - u) + (p/E')d (6) 

where E' is the effective Young's modulus taking into account the no-slip boundary condition on the upper and lower 
surface, which was measured to be E' — 4.2 MPa in a separate experiment where the rubber block was squeezed 
between two flat steel surfaces. Using (6) gives 

dp dp ds dp / Kd\ 
du ds du ds \ E' J 

or 

K* 
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where K* = dp/ds. Using (7) in Fig. [4] shows the normal contact stiffness as a function of the applied nominal 
contact pressure obtained from the measured p(s) relation with E' = 4.2 MPa (measured value) and E' = 4 MPa (to 
indicate the sensitivity of the result to E'). For very small contact pressures K* ps so that the denominator in (7) 
is « 1 (and K w K* as assumed in Ref. Q without proof) and the result is insensitive to E' as also seen in Fig. 
21 For large contact pressure the experimental data exhibits rather large noise (and great sensitivity to £"), which 
originates from the increasing importance of the compression of the rubber block for large contact pressure. That is, 
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FIG. 4: The normal contact stiffness as a function of the applied nominal contact pressure for a silicon rubber block (cylinder 
shape with diameter D = 3 cm and height d = 1 cm) squeezed against a road asphalt surface. The green and red lines are 
obtained from the measured p(s) relation using (7) with E' = 4.0 MPa and 4.2 MPa (see text) while the blue line is the theory 
prediction. 



for large pressures the denominator in (7) almost vanishes, which implies that a small uncertainty in the measured 
p(s) relation (which determines K*), or in E' , will result in a large uncertainty in K for large pressures. 

The blue curve in Fig. Uis the theory prediction which is obtained without any fitting parameter using the measured 
surface roughness power spectrum. For small contact pressure the contact stiffness obtained from the measured data 
is larger than predicted by the theory, but for nominal contact pressures typically involved in rubber applications 
(which are ~ 0.4 MPa as in tire applications, or higher in most other applications) the finite size effects are not 
important. 
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